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Abstract 

In this paper we announce a closed form analytical expression for the minimum attain- 
able attitude error (as well as the error rate) in a flexible beam by feedback control using 
co-located rate sensors. For simplicity, we consider a beam clamped at one end with an 
offset mass (antenna) at the other end where the controls and sensors are located. Both 
control moment generators and force actuators are provided. The results apply to any beam- 
like lattice-type truss, and provide the kind of performance criteria needed under CSI — 
Controls-Structures-Integrated optimization. 
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1. Introduction 

One of the challenges in the Design Challenge For Flexible Flight Structure Control 
System Design formulated in the inaugural paper on SCOLE [1] was to hold the antenna 
pointing error within ±0.02 degrees after slewing by appropriate feedback control. In this 
paper we derive a closed form expression for the minimal achievable mean square pointing 
error using co-located rate sensors. A slightly simplified form of the SCOLE article (which 
eliminates rigid-body modes) is used: a cantilevered beam with an offset mass where the 
controls — both c.m.g.'s and force actuators — and the rate sensors are located. Our results 
are in terms of continuum model parameters — the uniform Bernoulli version is used. The 
beam dynamics are given in Section 2. The main results are in Section 3. We note that a 
technique for deriving equivalent Bernoulli beam parameters for various types of trusses is 
described by Noor and Anderson in [4], Recently Noor and Russell [5] presented equivalent 
anisotropic Timoshenko beam models for beam-like lattice trusses with an arbitrary degree of 
modal coupling, which appear to yield excellent agreement with modal frequencies derived 
from finite element models. Our theory is able to handle these Timoshenko models, and 
moreover we can also use it for rigid-body modes, although they are not included here. Thus 
our results can be used for any beam-like lattice truss structure. 
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2. The Model 


We consider a uniform Bernoulli beam clamped at one end with an offset mass (antenna) 
at the other end which also houses the sensors and actuators. See Figure 1. We allow for both 
force actuators and moment actuators. The sensors are rate gyros. Because of the clamping at 
one end, no rigid-body modes are involved and hence no attitude sensors are needed. 

We allow bending in two mutually perpendicular planes containing the beam axis, as 
well as torsion in the plane perpendicular to the beam axis, all uncoupled. The continuum 
model (uniform Beroulli beam) dynamics can then be described by the following partial differ- 
ential equations (similar to those in [2, 3]). Let the beam extend along the z-axis, 0 < s < L, 
and let u^(s, t), u 6 (s, t), denote the bending displacements and u y ( j, r) the torsion angle 
about the beam axis. Let in the usual notation (cf. [1]), £/ e , EIq denote the flexural stiffness 
and G/ y the torsional rigidity. Let p denote the mass per unit area and A the cross-sectional 
area. Then we have: 


pA 
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<Pu^(s, 0 
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cp-u^s, f) 

" ari 
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+ ei 9 


d*u+(s, t) 

d? 
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-0, 0 < s < L \ 0 < t 

“0, 0<s<L;0<r 

0, 0<s<L;0<t 


with the clamped boundary conditions at s - 0 : 


V°>') = «e( 0 , t) » Ky( 0 .') - 0 

“*'(0, 0 - « e '(0, 0 - 0. 

The antenna center of gravity is located at 

(r x , r y ,L) . 

The distance from the beam tip to antenna center of gravity is denoted by 

M “ Vr* + . 
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Figure 1: Shuttle/ Antenna Configuration 
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The force balance equation at s - L yields 


u 6 (L, t) 

1 0 r x ♦ f x (t) EW(L,t) 

m u 0 (L, t) + - T w 

0 1 * f 2 (t) EI 9 Uq'(L, t) 

Uy(L,t) 

where m is the antenna mass and f\(), f 2 {) are the applied control forces. The torque 
balance equations yield 



EI+ 

AW 

0 = 

El 6 Uq(L, t) + / a tb + M(t) + r® 

flit) 


GlyUy(L,t) 

0 


uAL,t) + r x u {L, t) 
+ r 0 T 

ii e (L,f) + r y iiy(L, t) 


where the superdots indicate time derivatives and the primes the derivatives with respect to 
the spatial variable s ; ® denotes the vector cross-product and 0 ) the angular rate vector 

to = u^{L,t) , 

u V (L,0 

I a denotes the moment of inertia of the antenna about the beam tip (s - L) and finally, M(t) 
denotes the applied control moment. 

It is convenient to denote by 6(0 the boundary vector: 

u e (L, t) 

b(t) = u ;{L,t) . 

«e (L, t) 
iU t) 

The boundary rate vector would thus be 6(0 . Hence our sensor model is: 

v(r) = 6(0 + NJt) 

where we assume that N 0 (t) is white Gaussian noise with spectral density matrix d 0 I, where 
/ is the identity (5x5) matrix. Similarly we assume that the control actuators are also 
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characterized by additive white Gaussian noise. Denoting the applied control vector by u(f) . 

«2 (^. 0 

u(t) - u^{L,t) . 

«4 {L, 0 
us(L, t) 

we have 

fi(t) 

flit) - u(t) + N s (t) 

Mi t) 

where N s (t) is white Gaussian with spectral density d s . We shall also use M b to denote the 
actuator mass/inertia matrix 

m 0 0 0 mr x 

0 m 0 0 mr y 

M b - 0 0 

0 0 /„ 
mr x mr y 

where 

r\ ~r x r y 0 

K ” 7 « + ~ r x r y r x 0 

0 0 >x + T y 

where l a is the antenna moment of inertia about its center of gravity. For any control input 
u( ) (which must perforce be a “feedback” control, based on the sensor data v( )) the mean 

square pointing error is then expressed by: 

T t "I 

1 im ^ | f uAU tf dt + J u b (L, tf dt + |r| 2 J Uy(L, tf dt > 

T-*~ T l 0 0 o J 

and the mean square pointing rate is given by 
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. [ T T T 

1 im = \ J u^(L, tf dt + J u 0 (L, t) 2 dt + |r | 2 f i^(L, t f dt ’ . 
r ->~ U o o J 

From the results in [6] it follows that the minimal attainable mean square pointing error is 
given by 

a M~ b 'a* 

where 

a - row vector (1, 1, 0, 0, |r|) 
a * - transpose of a 


3. Main Results 

We need some notation first. The mean square attitude response, whatever the feed- 
back control used is defined by 

{ T T T 1 

J M'- *)’ dt + I «e('. *) 2 dt + |r| 2 f u (t, l ) 2 dt\. (3.1) 

o o o J 

This is recognized as the mean square displacement of the center of gravity of the antenna 

which is then also proportional to the mean square “pointing” error — see [1] for the 

relationships. 

Next let u denote any (vector) of control inputs — a constant “step” input: 


u 


“2 

“3 

«4 

«5 


(3.1) 


Solve the equations 


EW'(s) « 0 ’ 
EWXs) - 0 ' 

Gl w u^(s) - 0 . 


0 < s < L , 


(3.2) 
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subject to the end conditions 


El^u’;\L) - u, 
EI 6 uS\L) - u 2 
EI^u;(L) + u 3 - 0 
E! 6 u^L) + U4 - 0 
G/ v u^(L) + m 5 - 0 j 


(3.3) 


Note that the solution can be recognized as the steady-state response of the system to the 
step-input u, assuming that there is some damping. We only need to calculate die response to 
three specialized inputs: 

Calculate the response to u^(L) to the special case. Case 1, where: 

4, - 1 

4,-0, 2 £ 1 £ 5 . 

Calculate the response u e (L) to the special case. Case 2, where: 

4i - 0 

4 2 -l 

u 3 - 44 - 45 - 0 . 

Calculate the response Uy (L) to the special case, Case 3, where 

41 - 42 - 43 - 44 -0 

46 - 1 . 

Then the minimal achievable mean-square response whatever the choice of the feedback and 
whatever the mean-square control effort, is given by 

<d s d 0 (« 0 (L) 2 + Ufl(L) 2 + r 2 4y(L) 2 ) . (3.4) 

This is our main result. Unfortunately the derivation is beyond the scope of this report and 
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will be published elsewhere. To proceed further with (3.4) we calculate the solution of (3.2), 
(3.3) explicitly. Thus for any u, the solution is of the form 


where 


Thus for Case 1 we have 


and for Case 2 we have 


and for Case 3: 


«<,(*) - <*iS* + 


0 < s < L 

Uq(s) - + 

*2 J 2 . 

0 < s < L 

Uy(s) - C,J , 


0 < s < L 
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u # (L ) 2 - 


1? 


u 6 (L ) 2 


L? 

3 £/ e 


UyfL) 2 - 



Hence the mean-square attitude error 


VST UL + Jd- + 

* 0 13E/* + 3£/ 0 + G/ ¥ J • 


(3.5) 


Note the appearance of the noise parameters in (3.5) in product form. 

The technique for calculating the minimal mean square atttiude error in more complex 
models than that illustrated is the same: calculate the mean square step response (assuming 


457 



some damping) to unit step inputs. 

In conclusion we suggest this result (3.5) can be the basis for combined structures- 
controls optimization — CSI, since the required structural parameters can be calculated for a 
lattice truss from the material gage and physical dimensions as in [4, 5]. We omit the details 
of these calculations. 
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